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Highly  loaded  Low-Pressure  Turbine  (LPT)  blades  offer  opportunities  for  lower  costs 
and  higher  performance  but  suffer  from  unacceptable  endwall  losses.  This  has  motivated 
research  aimed  at  reducing  the  endwall  losses  by  manipulating  the  passage  vortex  which 
constitutes  the  primary  source  of  losses.  Passive  approaches  such  as  profile  and  endwall 
contouring  or  active  flow  control  did  reduce  endwall  losses  but  require  a  more  detailed 
understanding  of  the  underlying  flow  physics  to  be  better  optimized.  Endwall  research 
also  suffers  from  a  lack  of  temporal  and  spatial  resolution  of  the  existing  flow  data  which 
limits  the  depth  of  the  physical  analysis.  LPT  endwall  flow  experiments  at  the  Air  Force 
Research  Laboratory  revealed  the  mean  flow  features  but  did  not  resolve  the  unsteady 
flow  structures.  Therefore  the  decision  was  made  to  carry  out  endwall  simulations  with  a 
higher-order-accurate  Navier-Stokes  code  and  to  investigate  the  underlying  hydrodynamic 
instabilities.  In  this  paper  results  from  endwall  flow  simulation  for  the  L2F  geometry  with 
and  without  endwall  fillet  are  presented.  Both  a  turbulent  and  a  laminar  endwall  boundary 
layer  are  being  considered.  For  the  the  former  a  passage  vortex  is  generated  at  the  endwall. 

The  addition  of  the  fillet  suppresses  the  vortex.  For  the  laminar  endwall  boundary  layer, 
without  fillet  the  passage  vortex  is  missing.  This  interesting  result  indicates  a  strong 
dependence  of  the  flow  topology  on  the  endwall  boundary  layer  properties. 

I.  Introduction 

The  Low-Pressure  Turbine  (LPT)  produces  the  bulk  net  power  in  many  jet  engines.1  The  LPT-driven 
fan  is  responsible  for  up  to  80%  of  the  total  engine  thrust  and  the  LPT  can  constitute  as  much  as  one  third 
of  the  overall  engine  weight.2  Changes  in  LPT  efficiency  can  result  in  nearly  equal  changes  in  overall  engine 
efficiency.  Modern  LPTs  have  to  drive  larger  fans  at  lower  fan  speeds  and,  at  the  same  time,  have  to  be  less 
complex  and  lighter.  The  goal  in  modern  jet  engine  development  is  a  reduction  in  LPT  stage  solidity  without 
compromising  performance.  The  LPT  blade  count  is  limited  by  the  blade  loading.  A  higher  loading  would 
allow  for  a  reduction  of  the  overall  engine  cost  and  weight  and/or  a  performance  increase.  Ultimately,  the 
economics  of  air  travel  would  be  improved.  Towards  this  end  highly  loaded  LPT  blades  are  being  developed, 
with  loading  approaching  twice  the  currently  fielded  levels. 

LPTs  must  operate  efficiently  over  a  large  range  of  Reynolds  numbers.  Low  Reynolds  number  operating 
conditions  in  combination  with  aggressive  designs  can  lead  to  laminar  separation  from  the  suction  side,  which 
can  cause  significant  reductions  in  turbine  and  overall  engine  performance.  Highly  loaded  blades  are  typically 
front-loaded  to  counteract  laminar  separation  from  the  suction  side.  Front-loading  results  in  a  less  aggressive 
suction  side  pressure  recovery  and  thus  increases  stall  resistance  at  low  Reynolds  numbers.  Front-loading 
can  be  achieved  by  increasing  the  stagger  angle  (chord  angle  relative  to  the  incoming  flow)  or,  for  a  fixed 
stagger  angle,  by  thickening  the  leading  edge  (e.g.,  Korakianitis3  and  Korakianitis  and  Papagiannidis4). 
An  example  of  a  highly  loaded  design  is  the  L2F  profile  which  was  developed  at  the  Air  Force  Research 
Laboratory  (AFRL)  at  Wright- Patterson  Air  Force  Base. 
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Figure  1.  Sketch  of  endwall  flow  from  Langston.7 


Unfortunately,  for  the  aspect  ratios  encountered  in  typical  gas  turbine  engines,  front-loaded  airfoils  with 
high  stagger  angle  suffer  from  unacceptable  endwall  losses  (i.e.,  secondary  flow  losses).  For  that  reason, 
highly  loaded  profiles  which  promise  to  reduce  the  engine  parts  count  have  never  been  employed  in  production 
engines.  The  endwall  flow  is  highly  complex  (e.g.,  Sieverding,5  Wang  et  ah, 6  Langston7)  and  characterized 
by  the  passage  vortex  which  is  the  dominant  feature,  the  horseshoe  vortex,  the  counter  vortex,  and  the 
endwall  crossflow  among  others  (Fig.  1). 

Increased  endwall  losses  for  front-loaded  profiles  were  observed  by  a  large  number  of  researchers,  e.g., 
Weiss  and  Fottner,8  Zoric  et  ah, 9  and  Knezivici  et  al.10  For  example,  Praisner  et  al.11  found  high  endwall 
losses  for  the  Pack  DF  high  lift  design  which  has  a  25%  higher  loading  than  the  (conventional)  Pack  B 
blade.  The  front-loaded  profiles  in  these  earlier  investigations  were  designed  with  increased  stagger  angles 
compared  to  the  companion  aft-loaded  profiles.  Results  by  Lyall12  indicate  that  the  high  stagger  angle 
and  not  the  front-loading  is  responsible  for  the  endwall  losses.  High  stagger  profiles  tend  to  produce  a 
larger  “blockage”  for  the  incoming  endwall  boundary  layer  flow  (stronger  adverse  pressure  gradient)  thereby 
aggravating  endwall  separation  and  resulting  in  more  fluid  being  entrained  in  the  endwall  vortices. 

Endwall  losses  are  common  to  all  gas  turbines  stages  and  several  passive  flow  control  techniques  have  been 
proposed  to  alleviate  them  (for  a  review  see  Langston7),  among  them  boundary  layer  fences,  leading  edge 
bulbs  and  fillets  (e.g.,  Sauer  et  al.,13  Zess  and  Thole,14  Becz  et  al.,15  Saha  et  al.,16),  skewing  or  compound 
leaning  of  the  airfoil  (e.g.,  Harrison17),  and  endwall  contouring  (e.g.,  Rose,18  Harvey  et  al.,19  Hartland  et 
al.,20  Praisner21).  These  methods  aim  to  reduce  the  endwall  losses  by  manipulating  the  endwall  vortex 
structures  to  reduce  the  total  pressure  loss. 

Endwall  contouring  (or  profiling)  where  non-axisymmetric  endwalls  are  used  in  an  attempt  to  modify 
the  flow  has  been  investigated  extensively  for  high-pressure  turbines  (e.g.,  the  “Durham  cascade,”  Bagshaw 
et  al.22).  For  example,  Harvey  et  al.19  employed  non-axisymmetric  endwall  contouring  to  decrease  the 
cross-passage  pressure  gradient  which  weakened  the  passage  vortex  and  reduced  the  total  pressure  loss. 
McIntosh  et  al.23  employed  a  genetic  algorithm  for  optimizing  the  endwall  shape.  The  technology  has 
progressed  considerably  and  is  beginning  to  see  implementation  on  the  most  recent  engines  (e.g.,  Trent 
900).  Endwall  contouring  does,  however,  not  seem  to  work  sufficiently  for  highly  loaded  LPT  profiles.  For 
example,  Knezevici  et  al.10  and  Praisner  et  al.11  found  that  endwall  losses  remained  beyond  practical  limits 
for  the  high  lift  Pack  DF  profile,  even  after  implementing  non-axisymmetric  endwall  contouring.  Other 
approaches  such  as  profile  contouring  (or  combinations  of  different  approaches)  have  to  be  considered  to 
reduce  the  endwall  losses  of  high  lift  LPT  designs  to  acceptable  levels.  As  described  by  Langston,7  the 
concept  of  leading  edge  bulbs  is  to  strengthen  the  counter  vortex,  which  has  the  opposite  sense  of  rotation 
of  the  passage  vortex,  thus  weakening  the  overall  vortex  system  that  comprises  the  endwall  flow.  On  the 
other  hand,  Zess  and  Thole14  showed  that  a  leading  edge  fillet  can  reduce  or  eliminate  the  horseshoe  vortex 
that  forms  as  the  inlet  boundary  layer  separates.  Becz  et  al.15  accomplished  nearly  the  same  aerodynamic 
loss  reduction  with  both  a  leading  edge  bulb  and  a  fillet. 

The  AFRL  decided  to  investigate  profile  contouring  to  reduce  the  endwall  losses  (e.g.,  Lyall  et  al.24 
and  Sangston  et  al.25).  A  new  geometry  was  designed  that  blends  from  a  L2F  profile  inside  the  passage 
to  a  low-stagger  L2F-LS  airfoil  near  the  endwall.  The  “combined”  geometry  was  designated  L2F-EF.  This 
design  choice  was  made  based  on  the  observation  that  a  reduction  of  the  stagger  angle  near  the  endwall  would 
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weaken  the  end  wall  vortices  and  lower  the  endwall  losses.  Contouring  the  L2F  airfoil  at  the  endwall  to  obtain 
the  benefit  of  low  stagger  angle  generates  a  fillet  that  extends  out  from  the  pressure  surface,  primarily  in  the 
pitchwise  direction  toward  the  adjacent  suction  surface.  Experiments  for  Re=100,000  showed  an  approximate 
10%  reduction  of  the  total  pressure  loss  coefficient  for  the  contoured  L2F-EF  airfoil.12  Since  AFRL  is  carrying 
out  detailed  flow  field  measurements  for  this  configuration  it  was  chosen  for  the  present  simulations. 

The  profile  contouring  in  the  AFRL  design  did  reduce  endwall  losses  but  requires  a  more  detailed  un¬ 
derstanding  of  the  endwall  flow  physics  to  be  better  optimized.  Reynolds- Averaged  Navier-Stokes  (RANS) 
calculations  that  were  carried  out  during  the  design  process24  provided  reasonably  accurate  mean  flow  predic¬ 
tions  but  suffered  from  the  following  shortcomings:  By  design  most  RANS  models  cannot  capture  transition 
which  has  been  found  to  play  a  dominant  role  in  LPT  flows.  By  assuming  fully  turbulent  flow,  the  laminar 
separation  bubble  that  can  naturally  occur  on  the  suction  surface  is  suppressed.  Laminar  separation  bub¬ 
bles  not  only  transition  the  flow.  They  also  act  as  natural  “oscillators”  and  shed  coherent  structures  which 
introduce  a  large  amount  of  unsteadiness  into  the  flow.  In  RANS  calculations  this  unsteadiness  is  inhibited 
because  the  bubble  is  not  captured.  In  addition,  since  RANS  models  are  derived  under  the  assumption 
that  all  unsteady  motion  is  modeled,  they  typically  introduce  a  large  amount  of  turbulent  eddy  viscosity 
which  dampens  out  all  or  most  of  the  unsteadiness.  As  a  result  not  all  of  the  relevant  flow  physics  are 
captured  in  RANS  calculations.  Simulations  of  high-lift  LPT  endwall  flows  at  low  Reynolds  numbers  are 
quite  challenging  primarily  due  to  the  intricate  interaction  of  transition  and  separation.  Direct  Numerical 
Simulations  (DNS;  the  highest  fidelity  approach)  or  Implicit  Large  Eddy  Simulations  (ILES)  that  capture 
the  essential  flow  physics  at  reduced  computational  cost  appear  to  be  the  most  feasible  alternatives.  As  far 
as  LPT  simulations  with  endwall  are  concerned,  P.  Tucker  and  co-workers  at  Cambridge  University  (UK)  in 
collaboration  with  Rolls-Royce  appear  to  be  the  only  ones  that  have  carried  out  a  time-accurate  simulation  of 
a  LPT  cascade  with  endwall.  Some  flow  visualizations  are  available  online  but  to  the  authors  best  knowledge 
none  of  the  results  have  been  published. 

In  this  paper  results  from  time-resolved  high-resolution  endwall  flow  simulation  for  the  L2F  and  L2F-EF 
geometries  are  presented.  First  the  methods  employed  for  generating  the  computational  grid  and  the  tur¬ 
bulent  endwall  boundary  layer  are  explained.  Then  simulation  results  for  a  turbulent  and  a  laminar  endwall 
boundary  layer  are  presented  and  the  total  pressure  loss  coefficients  for  the  various  cases  are  compared.  The 
paper  concludes  with  a  brief  summary  and  conclusions. 

II.  Methodology 


A.  Simulated  Cases 

Simulations  were  carried  out  for  the  L2F  blade  which  was  developed  at  AFRL  for  investigating  the  low- 
Reynolds  number  and  high  lift  LPT  aerodynamics.26  The  design  intent  was  to  have  good  low- Reynolds 
number  performance  compared  to  the  Pratt  and  Whitney  Pack  B  blade  and  high  aerodynamic  lift.  The 
axial  blade  chord  and  span  in  the  AFRL  experiments25  were  Cx=6 in  and  H=3.5CX  (21in),  respectively.  The 
pitchwise  spacing  was  1.221  Cx  and  thus  38%  larger  than  for  the  Pack  B  blade.  The  L2F  blade  has  a  design 
inflow  angle  of  35deg  and  an  exit  angle  of  -58.12deg.  The  stagger  angle  is  34.60deg.  In  addition,  the  L2F-EF 
geometry  (L2F  with  endwall  fillet)  was  considered  which  blends  from  a  L2F  airfoil  at  mid- span  to  a  L2F-LS 
at  the  endwall.24  The  blending  extends  over  9.5%  of  the  total  span  (=0.3325 Cx)  at  both  endwalls.  The 
simulations  were  set  up  according  to  the  AFRL  experiments  where  upstream  wakes  and  surface  roughness 
were  not  considered.  This  deliberate  simplification,  when  compared  to  the  environment  in  real  jet  engines, 
allows  for  an  easier  extraction  of  the  flow  physics  and  a  comparison  with  the  experimental  data.  Only  half 
of  the  blade  span  was  simulated  and  symmetry  conditions  were  enforced  at  mid-span.  The  Reynolds  number 
based  on  axial  chord,  Cx ,  for  the  simulations  was  Re=100,000  and  thus  the  same  as  in  the  experiments  and 
earlier  RANS  calculations.24,25 

B.  Endwall  Boundary  Layer  State 

The  endwall  in  the  experiment  had  an  elliptic  leading  edge  with  (in  the  axial  direction)  major  semiaxis  of 
0.25 Cx  (1.5in)  and  (in  the  wall-normal  direction)  minor  semiaxis  of  0.0417 Cx  (0.25in).  Measured  in  the  axial 
direction,  the  flat  plate  leading  edge  was  located  at  x=-3.95 8CX  (23.75in)  where  the  minus  sign  indicates 
a  location  upstream  of  the  inflow  plane  of  the  cascade  (x=0).  The  corresponding  streamwise  distance  is 
s=-3. 958/cos  35deg=-4. S33CX  (29in).  The  displacement  thickness,  5*  =  f(  1  —  v/v^ )dy,  for  a  laminar  flat 
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plate  boundary  layer  is  5*  =  1.721  s/Re®'5.  For  s  =  4.833 Cx,  Res *  =  1.721  y/s/Cx  Re  =  1,200.  At  this 
Reynolds  number  (for  the  experimental  freestream  turbulence  intensity,  FSTI=1%)  the  boundary  layer  (at 
the  cascade  inlet  plane)  can  be  expected  to  be  turbulent. 

In  the  RANS  calculations  by  Lyall  et  al.24  a  turbulent  boundary  layer  with  a  thickness  of  £99=0.02577=0.087567 
was  prescribed  at  the  inflow  boundary  which  was  located  at  x=-  1ACX.  The  following  approximations  hold 
for  a  turbulent  flat  plate  boundary  layer: 

e  0.37s  _  0.046s  ^  0.036s  1  0.0592 

SM  =  Tte F  •  s  “flJF  '  ,=  &F  •and  Rip-  (1) 

Accordingly,  the  distance  to  the  “virtual  leading  edge”  (measured  in  the  streamwise  direction)  for  the  Lyall 
et  al.  calculations  is  s/Cx=(l/0.37  Sqq/Cx  77e°-2)1/°-8=2.933.  Relative  to  the  inflow  plane  of  the  cascade,  the 
virtual  leading  edge  is  located  at  £=-2.933xcos  35deg-1.4=-3. 802(7^.  From  the  inflow  velocity  profiles  of  the 
Lyall  et  al.  calculations24  the  displacement  and  momentum  thickness  were  computed  as  £*=0.01254367  and 
0=0.00923536^  which  corresponds  to  5/6^=3.504  (x =-4.2 70CX)  and  s/Cx=3.247Cx  (£=-4.0606^). 

The  endwall  in  the  experiments  is  segmented  (i.e.,  modular)  and  different  upstream  lengths  are  possible. 

The  shortest  possible  upstream  extent  of  the  endwall  is  0.683C7  (4.1in)  in  the  axial  direction  and  0.834(7^ 

(5in)  in  the  streamwise  direction.  The  corresponding  Reynolds  number  based  on  displacement  thickness  is 
Res*  =497.  At  this  Reynolds  number  the  boundary  layer  (at  the  cascade  inlet  plane)  is  likely  laminar. 

C.  Navier-Stokes  Code 

A  research  computational  fluid  dynamics  code  that  was  developed  in  our  laboratory  was  employed  for  the 
present  investigations.27  The  code  solves  the  compressible  Navier-Stokes  equations  in  the  finite  volume 
formulation.  The  convective  terms  were  discretized  with  a  ninth-order-accurate  upwind  scheme  (unless 
stated  otherwise).  The  viscous  terms  were  discretized  with  a  fourth-order-accurate  scheme.  A  second-order- 
accurate  implicit  Adams-Moulton  method  was  employed  for  time  integration. 

D.  Normalization 

Length  scales  were  made  dimensionless  with  the  axial  chord  length,  Cx.  Velocities  were  made  dimensionless 
with  the  cascade  inflow  velocity,  v^.  Density  was  made  dimensionless  with  the  inflow  density,  poo,  and 
pressure  was  made  dimensionless  with  Poovlo-  Time  was  made  dimensionless  with  Cx/voq.  The  reference 
Mach  number  was  M= 0.1.  This  Mach  number  is  small  enough  to  satisfy  the  incompressible  flow  assumption 
without  negatively  affecting  the  convergence  characteristics  of  the  implicit  time-integration  scheme.  The 
Prandtl  number  was  Pr =0.72. 

E.  Computational  Grids 

Two  different  computational  grids  were  generated  with  a  Poisson  grid  generator28  (Fig.  2).  The  block 
boundaries  are  outlined  in  red.  Both  grids  are  periodic  in  the  pitchwise  direction.  The  virtual  flat  plate 
leading  edge  is  indicated  by  a  bold  vertical  line.  The  first  grid  (grid  1)  consists  of  four  blocks  and  extends 
upstream  to  £=-2.9536.  This  grid  was  created  for  the  case  where  the  endwall  leading  edge  is  at  £=-3.958 
and  the  endwall  boundary  layer  is  turbulent.  The  second  grid  (grid  2)  consists  of  two  blocks  and  extents 
upstream  to  £=-0.5.  This  grid  was  designed  for  the  case  where  the  endwall  leading  edge  is  at  £=-0.6833  and 
the  endwall  boundary  layer  is  laminar.  For  both  grids,  streamwise  grid  line  stretching  was  employed  near 
the  outflow  to  dissipate  wake  structures.  Also  for  both  grids,  singular  points  are  obtained  where  the  central 
O-grids  connect  with  each  other.  The  block  grid  resolutions  are  provided  in  Tab.  1.  For  block  1  the  number 
of  cells  in  the  wall- normal  direction  is  larger  for  grid  1  compared  to  grid  2.  This  was  mainly  done  to  obtain 
a  reasonable  pitchwise  grid  resolution  for  block  4  with  a  cell  aspect  ratio  (streamwise  to  pitchwise)  close  to 
one. 

For  the  3-D  simulations  (without  endwall  fillet)  the  2-D  grids  were  “extruded”  in  the  spanwise  direction. 

For  both  cases  K= 256  cells  were  employed  in  the  spanwise  direction.  The  total  number  of  cells  was  15.2 
million  for  the  turbulent  case  (grid  1)  and  9.27  million  for  the  laminar  case  (grid  2).  For  comparison,  713 
thousand  cells  were  employed  for  the  calculations  by  Lyall  et  al.24  The  wall- normal  grid  line  spacing  at  the 
endwall  was  Az  =  lO-4^.  A  grid  line  distribution  that  provides  grid  line  clustering  near  the  endwall 

^o  =  0  (2) 
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Figure  2.  Details  of  computational  grid  for  turbulent  case  (left,  grid  1)  and  laminar  case  (right,  grid  2). 


Block 

Turbulent  case,  grid  1 

Laminar  case,  grid  2 

1 

400x120 

400x80 

2 

60x100 

3 

20x100 

30x140 

4 

160x20 

total 

59,200 

36,200 

Table  1.  Number  of  cells  per  block. 


Zl 


Zk 


A  * 


Zk- 1  + 


0.1 


(k-K 

V2-K 


+  1 


(Zfc-1  -  Zk- 2) 


for  k  >  2 


(3) 

(4) 


and  blends  towards  an  equidistant  spacing  near  mid-span  was  employed  in  the  spanwise  direction.  The  near- 


Figure  3.  Near- wall  grid  resolution  in  wall  units  (From  left  to  right:  Wall- normal/  circumferential/  spanwise- 
radial). 


wall  grid  resolution  for  the  coarser  grid  (grid  2)  is  provided  in  Fig.  3.  The  wall- normal  near- wall  grid  line 
spacing  is  less  than  one.  The  circumferential  (streamwise)  grid  line  spacing  in  wall  units  at  the  wall  is  less 
than  80  and  the  spanwise  grid  line  spacing  is  about  100.  Georgiadis  et  al.29  recommend  a  grid  resolution  of 
50  <  Ax+  <  150,  A y+  <  1,  15  <  Az+  <  40  for  Large  Eddy  Simulations  (LES).  Accordingly  the  streamwise 
and  wall- normal  resolution  of  the  present  grid  is  sufficient.  To  meet  the  Az+  requirement,  the  spanwise  grid 
resolution  would  have  to  be  increased  2.5  times.  The  near- wall  resolution  for  the  grid  1  (turbulent  case)  is 
slightly  higher  and  not  shown. 

For  the  L2F-EF  geometry  the  LPT  coordinates  of  the  computational  grid  by  Lyall  et  al.24  were  splined 
in  the  ^-direction  and  z=const.  slices  were  extracted  (Fig.  4).  For  each  slice  a  2-D  grid  was  generated  with 
the  Poisson  grid  generator.  In  between  the  z-slices  shown  in  Fig.  4  the  grid  points  were  interpolated  (in  the 
spanwise  direction)  using  either  Hermite  or  Bezier  splines.  The  Bezier  splines  did  (by  design)  not  exactly 
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Figure  4.  Slices  of  constant  z  for  L2F-EF  geometry. 


pass  through  the  control  points  but  resulted  in  an  overall  higher  grid  quality  (Fig.  5).  Considering  the 
relatively  low  geometry  information  content  of  the  source  grid  (limited  number  of  z-planes)  the  fact  that  the 
Bezier-spline-based  grid  missed  the  control  points  by  a  very  small  amount  (this  was  confirmed  visually)  was 
determined  acceptable.  In  the  following,  the  grid  with  fillet  derived  from  grid  1  (turbulent  endwall  boundary 
layer)  is  referred  to  as  grid  3  and  the  grid  derived  from  grid  2  (laminar  endwall  boundary  layer)  is  referred 
to  as  a  grid  4.  The  number  of  cells  per  block  remains  unchanged. 


Figure  5.  From  left  to  right:  1)  Lyall  et  al.24  grid;  2)  Hermite  spline  grid;  3)  Bezier  spline  grid. 


F.  Boundary  Conditions 

Walls  were  considered  as  adiabatic.  Non-reflecting  boundary  conditions30  were  applied  at  the  inflow  and 
outflow  boundary.  Symmetry  was  enforced  at  mid-span  (only  half  of  the  passage  span  was  simulated). 
Periodicity  conditions  were  applied  in  the  pitchwise  direction. 

1.  Turbulent  Endwall  Boundary  Layer 

The  purpose  of  block  4  for  grid  1  (Fig.  2)  was  to  obtain  an  equilibrium  turbulent  boundary  layer  at  the 
inflow  to  block  3.  A  Blasius  velocity  profile  with  Res  =  264,  810  was  prescribed  at  the  inflow  boundary  of 
block  4.  A  forcing  slot  of  length  Lx  was  added  near  the  inflow  boundary  of  block  4  to  introduce  disturbances 
into  the  boundary  layer  (the  forcing  slot  is  highlighted  in  black  in  Fig.  6).  Wall- normal  blowing  and  suction 
with 

w  forcing  =  ax^  Ai  cos  (atx  +  /%  -  UJit)  (5) 

i 

and  axial  and  pitchwise  wavenumbers  of  —  2tt/X ^x  and  =  27r/A was  applied  for  -2.8901<  x  <-2.7631. 
A  half-sine  shaped  amplitude  distribution  was  prescribed  in  the  axial  direction,  ax  =  sin((x -f  2.8901  )/Lx  tt). 
It  was  decided  to  force  a  Tollmien-Schlichting  wave  and  a  steady  3-D  disturbance.  Diagrams  found  in  Baines 
et  al.31  provided  guidance  for  the  selection  of  the  period  and  wavelength  of  the  traveling  wave.  Several 
parameter  combinations  were  tried.  The  combination  in  Tab.  2  resulted  in  a  rapid  breakdown  of  the  endwall 
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i 

A 

^ X 

Xy 

U) 

A 

1 

2 

0.2 

0.1221  x  cos35deg 

A 

A 

2ir 

0.1 

0.1 

cos  35 deg 
-A 

sin  35 deg 

A 

0.1 

0 

sin  35deq 

cos  35 deg 

Table  2.  Forcing  parameters. 


boundary  layer.  A  visual  representation  of  the  endwall  boundary  layer  is  provided  in  Fig.  6.  For  the  ninth- 
order-accurate  discretization  a  larger  wavenumber  interval  is  resolved  (compared  to  the  fifth-order-accurate 
discretization).  Time-  and  pitchwise-averaged  velocity  profiles  near  the  outflow  boundary  of  block  4  are 


Figure  6.  Block  4  endwall  boundary  layer.  Iso-surfaces  of  Q= 1  colored  by  velocity  magnitude  obtained  with 
fifth-  and  ninth-order-accurate  scheme  (left  and  right). 

shown  in  Fig.  7.  As  a  reference,  the  relationships  for  the  viscous  sublayer,  U+  =  ?/+,  and  for  the  log-layer, 


05-scheme  09-scheme 


Figure  7.  Block  4  velocity  profiles. 

u+  =  5  +  1/0.41  In y+  were  included.  For  x=-1.4,  Res= 312,000  and  c/=0. 00471  are  obtained  and  the  profiles 
are  scaled  according  to  u/u^  =  yj  Cf  j 2  x  u+  and  z/Cx  =  y/2/cf/Re  x  z+.  In  the  present  ILES  only 
the  larger  flow  structures  are  resolved  and  numerical  diffusion  dampens  out  much  of  the  high  wavenumber 
spectral  content.  As  a  result  the  near- wall  Reynolds- stresses  are  underpredicted.  The  skin- friction  coefficient 
is  shown  in  Fig.  8.  Also  included  are  reference  lines  for  a  laminar  ( Cf  =  (Cx/s)°'50. 664/Re0  5)  and  a  turbulent 
boundary  layer  (cf  =  (Cx/s)°-20. 0592/ Re0  2)  starting  from  the  endwall  leading  edge.  Downstream  of  the 
forcing  slot  the  skin-friction  coefficient  displays  a  mild  oscillation  which  can  be  attributed  to  the  forced 
steady  3-D  mode  which  is  decaying  in  the  streamwise  direction.  Overall  the  skin  friction  coefficient  is  too 
low  although  the  mismatch  is  less  severe  for  the  more  accurate  discretization. 

The  displacement  and  momentum  thickness  were  computed  for  four  downstream  locations  (Tab.  3).  Also 
included  in  the  table  is  the  shape  factor  H  =  5* /0.  As  a  reference,  from  the  Lyall  et  al.24  data  =0.012, 
0=0.0091,  and  77=1.36  were  extracted  for  x=-lA.  The  relationships  for  the  turbulent  boundary  layer  (Eqs. 
1)  provide  5* =0.011,  0=0.0089,  and  77=1.27.  Compared  to  these  reference  data,  both  displacement  and 
momentum  thickness  for  the  present  simulations  are  too  too  small.  The  difference  is  less  for  the  ninth-order 
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Figure  8.  Block  4  skin-friction  coefficient  for  Re  =  100,  000. 


accurate  scheme.  The  shape  factor  is  smaller  than  2.59  (Blasius  boundary  layer)  but  still  larger  than  what 
would  be  expected  for  a  turbulent  boundary  layer  ( H  «  1.3).  Although  the  turbulent  boundary  layer  in  the 


x/Cx 

S*/Cx 

e/cx 

H 

fifth-order-accurate  scheme 

-1.303 

0.0104 

0.00659 

1.57 

-0.922 

0.0109 

0.00687 

1.58 

ninth-order-accurate  scheme 

-1.303 

0.0103 

0.00690 

1.50 

-0.922 

0.0112 

0.00764 

1.47 

Table  3.  Block  4  boundary  layer  properties. 

present  simulations  does  not  accurately  match  the  reference  data  the  approach  appears  justified  because  of 
the  considerable  compute  time  savings  compared  to  a  fully-resolved  DNS. 

2.  Laminar  Endwall  Boundary  Layer 

The  smallest  possible  upstream  extent  of  the  endwall  from  the  cascade  inflow  in  the  experiment  is  x=-0.683. 
The  inflow  boundary  for  grid  2  is  at  x=-0.5.  The  difference  is  0.183  (in  the  axial  direction)  or  8=0. 183/cos 
35deg=0.223  in  the  streamwise  direction.  Accordingly,  a  Blasius  velocity  profile  with  Res= 22,300  was 
prescribed  at  the  inflow  boundary  of  block  1. 


G.  Total  Pressure  Loss  Coefficient 

With  the  chosen  non-dimensionalization  the  reference  (inlet)  dynamic  and  total  pressure  are 


Qref  =  \ Pref  («re/  +  vref  +  wlef )  =  \ 


P0,ref  —  Pref  T  Qref  — 


The  inflow  and  outflow  dynamic  and  total  pressure  are 


1 


7M2 


—  —  (  2  -i-  2  -i-  2  ) 

Qin,out  2  Pin, out  \^in,out  '  ^ in, out  '  ^ in, out) 

^  ^  ^  _  Pin,out^-'in,out  {  ^ 

P0,in,out  —  Pin, out  T  Q_in,out  —  yTf 2  '  Qin,out  • 
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The  total  pressure  coefficient  is  defined  as 


PO  P0,ref 

Qref 


The  mass-averaged  total  pressure  coefficient  was  approximated  as 


=  ff  p\v\cpto  dydz  _  E»  P\v\cp,ovi 
P’°  ff  P\v\dydz  ~  Y,iP\v\Vi 


(10) 


(11) 


where  the  sums  were  taken  over  the  cell  volumes,  V^,  directly  adjacent  to  the  inflow/outflow  boundary.  The 
same  procedure  was  applied  for  mass-averaging  the  dynamic  pressure.  The  total  pressure  loss  coefficient  was 
computed  as25 


Y, 


tot  ^ p,0,in  ^p,0,out  • 


(12) 


H.  Location  of  Slices  for  Comparison  with  Sangston  et  al.25  Measurements 

Sangston  et  al.25  acquired  Particle  Image  Velocimetry  (PIV)  data  in  an  outflow  plane  at  x=l.bCx  and  in 
four  planes  normal  to  the  chord  line  (Fig.  9).  The  L2F  stagger  angle  is  A=34.60deg.  The  origin  of  the 


Figure  9.  Slices  where  PIV  data  was  obtained  by  Sangston  et  al.25  Fig.  4a  from  Sangston  et  al.,25  coordinate 
systems,  and  3-D  representation. 


chord-normal  slices  is  at  xo=0.25,  0.5,  0.75,  1  and  yQ  =  (1  —  x0)tanA.  The  surface  normal  of  the  planes  is 
n  =  [cos  A,  —  sinA,0].  The  coordinate  transformation  from  (x,  y ,  z)  to  {xf,  y' ,  z')  is 


x'  =  (x  —  x0)  cos  A  —  (y  —  yQ)  sin  A  , 
y'  =  {x-  x0)  sin  A  +  {y  -  yQ)  cos  A  , 


and  z’—z.  The  vorticity  in  the  x'  direction  is 


uxf  =  uox  cos  A  —  cjy  sin  A  . 


(13) 

(14) 

(15) 


(16) 


III.  Results 

A.  Turbulent  Endwall  Boundary  Layer 

1.  Straight  Blade  (No  Fillet) 

Iso-surfaces  of  the  Q  vortex  identification  criterion  at  f=6  in  Fig.  10  illustrate  how  the  turbulent  endwall 
boundary  layer  interacts  with  the  blade.  The  vortex  identification  criterion,32 

Q  =  \  {WijWij  -  SijSij)  ,  (17) 

indicates  areas  where  rotation  dominates  strain.  The  endwall  region  near  the  suction  side  corner  is  relatively 
“quiet”  while  a  large  amount  of  flow  structures  is  seen  near  the  pressure  side  corner.  The  boundary  layer  on 
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Figure  10.  Instantaneous  iso-surfaces  of  Q  —  10. 


the  LPT  blade  (away  from  the  wall)  is  laminar  from  the  leading  edge  and  remains  laminar  on  the  pressure 
side.  The  suction  side  boundary  layer  separates  laminar  approximately  at  the  beginning  of  the  “uncovered” 
turning.  The  separated  boundary  layer  transitions  rapidly  and  reattaches  to  the  blade.  Also  noticeable  is 
the  pitchwise  periodicity  of  the  endwall  structures  at  the  inflow  of  block  2  that  results  from  the  relatively 
small  pitchwise  extent  of  block  4  (one  fifth  of  the  blade  spacing). 

In  Fig.  11  skin- friction  lines  and  streamlines  computed  from  the  time-average  are  overlaid  with  an  instan¬ 
taneous  flow  visualization.  The  skin-friction  lines  reveal  a  saddle  point  upstream  of  the  leading  edge  which 


Figure  11. 


Instantaneous  iso-surfaces  of  Q  =  10,  skin- friction  lines  (red  lines),  and  streamlines  (green  lines). 
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indicates  a  three-dimensional  (3-D)  separation.  In  between  the  blades  a  line  of  separation  is  seen  which  must 
be  associated  with  the  passage  vortex.  The  streamlines  (computed  from  a  z  =  0.05  slice  of  the  time-averaged 
data)  are  strongly  curved  at  the  beginning  of  the  passage.  According  to  the  radial  momentum  equation , 


dp  pv2 
dr  r 


(18) 


the  curvature  of  the  inviscid  flow  in  the  cascade  determines  the  radial  pressure  gradient  (between  the  pressure 
and  suction  side).  The  pressure  gradient  results  in  a  strong  endwall  crossflow  from  the  pressure  side  to  the 
suction  side  (Fig.  11). 


Figure  12.  Time-average.  Iso-surfaces  of  Q  =  10. 

Visualizations  of  the  time- averaged  flow- field  (all  time-averages  were  obtained  for  6  <  t  <  12)  are  provided 
in  Fig.  12  -  14.  The  passage  vortex,  which  is  hard  to  make  out  in  the  instantaneous  flow  visualizations,  clearly 
shows  up  in  the  time-averaged  data.  In  addition  to  the  passage  vortex,  a  number  of  streamwise  vortices  of 
varying  intensity  are  originating  from  the  endwall  (Fig.  12).  The  strong  pitchwise  endwall  flow  may  lead  to 
the  downstream  growth  of  steady  crossflow  vortices  (in  some  instances  the  skin-friction  lines  are  almost  at  a 
90deg  angle  with  respect  to  the  streamlines;  Fig.  11).  The  crossflow  vortices  may  be  “seeded”  by  numerical 
noise  or  by  numerical  inaccuracies  associated  with  the  singular  grid  point.  To  address  the  latter  concern, 
the  singular  point  could  be  moved  farther  upstream  where  the  wall  skin-friction  lines  follow  the  inviscid 
streamlines  or  a  different  grid  topology  could  be  conceived  where  the  singular  points  are  placed  at  entirely 
different  locations. 


The  dominant  structures  on  the  suction  side  are  a  vortex  associated  with  the  3-D  trailing  edge  separation 
near  the  blade  root  (which  results  in  a  line  of  separation  of  the  skin-friction  lines)  and  a  very  narrow  corner 
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vortex  (Fig.  13).  On  the  pressure  side  a  short  corner- vortex  can  be  associated  with  a  complicated  3-D  flow 
separation  near  the  leading  edge.  A  comparison  of  the  passage  vortex  location  with  the  PIV  measurements 


Figure  14.  Location  of  Passage  Vortex  (PV)  and  Suction  Side  Horseshoe  Vortex  (SHV).  Fig.  9a  from  Sangston 
et  al.25  overlaid  with  Q= 10  iso-contours  from  present  simulation. 


by  Sangston  et  al.25  is  provided  in  Fig.  14.  Compared  to  the  experiment  the  passage  vortex  is  slightly  closer 
to  the  suction  side  of  the  blade.  The  present  simulations  show  no  unique  suction  side  vortex  and,  therefore, 
a  comparison  with  the  suction  side  horseshoe  vortex  in  the  experiment  is  not  possible. 


1 
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Figure  15.  Wall  pressure  coefficient  iso-contours  and  skin-friction  lines. 


Figure  15  provides  wall-pressure  iso-contours  and  skin- friction  lines.  The  wall  pressure  increases  rapidly 
at  the  end  of  the  laminar  separation  bubble  on  the  suction  side  of  the  blade.  Near  the  leading  edge,  the 
corner  vortices  raise  the  pressure  on  the  suction  side  and  lower  the  pressure  on  the  pressure  side.  In  Fig.  16 
iso-contours  of  the  vorticity  in  the  direction  of  the  cascade  outflow  angle, 


uoout  =  cos  60 deg  —  uoy  sin  60 deg  , 

are  shown.  The  vorticity  vector  was  computed  as 


^ X 

d/dx 

u 

Ldy 

= 

d/dy 

X 

V 

_ 

d/dz 

w 

The  passage  vortex  is  stronger  and  more  coherent  than  the  vortices  on  the  suction  side  of  the  blade. 


(19) 


(20) 
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Figure  16.  Iso-contours  of  vorticity  magnitude,  uout‘ 


Figure  IT.  Vorticity  iso-contours,  ujxt,  and  in-plane  velocity  vectors.  Top:  PIV  measurements  by  Sangston  et 
al.25  Bottom:  Present  results. 

Iso-contours  of  the  chord- wise  vorticity,  av,  for  the  four  locations  indicated  in  Fig.  9  are  provided  in  Fig. 
17.  The  agreement  with  the  the  PIV  measurements  by  Sangston  et  al.25  is  adequate. 


Figure  18.  Iso-contours  of  ujout  at  x=1.5Cx •  PIV  measurements  by  Sangston  et  al.25  and  present  simulation. 


Iso-contours  of  the  streamwise  vorticity,  wout,  at  x=1.5Cx  in  Fig.  18  show  adequate  agreement  with  the 
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triple  probe  data  by  Sangston  et  al.25  The  present  near  wall  results  are  not  in  agreement  with  the  PIV 
measurements  by  Sangston  et  al.25  This  may  again  be  attributed  to  a  lack  of  grid  resolution  (the  grid 
resolution  is  degrading  quickly  downstream  of  the  blade,  see  Fig.  2). 

2.  Blade  With  Fillet 


Figure  19.  Instantaneous  iso-surfaces  of  Q  —  10. 


Next  the  AFRL  LPT  geometry  with  endwall  fillet  (L2F-EF)  was  considered.  Instantaneous  flow  visu¬ 
alizations  are  provided  in  Fig.  19.  The  endwall  region  near  the  pressure  side  of  the  blade  is  populated  by 
elongated  streamwise  structures  that  are  aligned  parallel  to  each  other  (again  indicating  an  instability).  The 
skin-friction  lines  between  the  blades  (Fig.  20)  do  not  feature  a  line  of  separation  and  the  z  =  0.05  streamlines 
follow  the  suction  side  contour  better  than  for  the  straight  blade  (Fig.  11). 
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Figure  20.  Instantaneous  iso-surfaces  of  Q  =  10.  Red  lines:  Skin-friction  lines  for  time  average.  Green  lines: 
Streamlines  at  z  =  0.05  for  time-average. 


Figure  21.  Time-average.  Iso-surfaces  of  Q  =  10. 

Visualizations  of  the  time- averaged  flow  field  in  Fig.  21  reveal  that  the  passage  vortex  is  significantly 
weaker  than  for  the  straight  blade  (Fig.  12).  This  is  in  agreement  with  the  measurements.25  The  passage 
vortex  trajectory  does,  however,  not  match  up  with  the  trajectory  obtained  from  the  PIV  measurements25 
(Fig.  22).  Skin- friction  lines  for  the  pressure  surface  of  the  fillet  (Fig.  23)  indicate  attached  flow  while  the 
skin- friction  lines  for  the  straight  blade  (Fig.  16)  indicate  separation  and  reverse  flow.  Iso-contours  of  the 
streamwise  vorticity  (Fig.  24)  provide  additional  confirmation  of  a  weakened  passage  vortex  compared  to  the 
straight  blade  (Fig.  16). 
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Figure  22.  Location  of  Passage  Vortex  (PV)  and  Suction  Side  Horseshoe  Vortex  (SHV).  Fig.  9a  from  Sangston 
et  al.25  overlaid  with  Q=10  iso-contours  from  present  simulation. 
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A  comparison  with  PIV  measurements  by  Sangston  et  al.25  is  provided  in  Fig.  25.  The  data  are  in 
qualitative  agreement.  In  the  experiment  the  suction  side  geometry  of  the  L2F-EF  was  changed  near  the 
endwall:  For  ease  of  manufacturing,  the  fillet  was  machined  as  an  extra  piece  that  was  added  to  the  L2F 
geometry.  As  a  result  the  recession  on  the  suction  side  of  the  L2F-EF  near  the  endwall  was  omitted  in  the 
experiment  which  may  have  a  slight  effect  on  the  flow  field.  Vorticity  iso-contours  at  the  outflow  boundary 
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Figure  25.  Vorticity  iso-contours,  cjx/,  and  in-plane  velocity  vectors.  Top:  PIV  measurements  by  Sangston  et 
al.25  Bottom:  Present  results. 

are  shown  in  Fig.  26. 
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Figure  26.  Iso-contours  of  uo0ut  at  x=1.5Cx 


3.  Comparison 

Direct  comparisons  of  the  results  for  the  L2F  and  L2F-EF  geometry  are  provided  in  Figs.  27  &  28.  The 
following  observations  can  be  made  regarding  the  endwall  (Fig.  27):  With  fillet  the  line  of  separation  (as 
seen  in  the  skin-friction  lines  for  the  straight  blade)  and  the  pressure  dip  (pressure  iso-contours)  associated 
with  the  passage  vortex  disappears  and  the  streamlines  at  z= 0.05  are  better  aligned  in  the  outflow  direction. 
Considering  the  blade  suction  surface,  with  fillet  the  spanwise  width  of  the  trailing  edge  separation  at  the 
endwall  is  reduced  (Fig.  28). 
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Figure  27.  Skin-friction  lines,  wall-pressure  iso-contours,  and  streamlines  at  z= 0.05.  Red  lines:  Straight  blade. 
Green  lines:  Blade  with  fillet. 


Figure  28.  Suction  side  skin-friction  lines  and  wall-pressure  iso-contours.  Red  lines:  Straight  blade.  Green 
lines:  Blade  with  fillet. 


B.  Laminar  End  wall  Boundary  Layer 

1.  Straight  Blade  (No  Fillet) 

When  the  approach  flow  endwall  boundary  layer  is  laminar,  traveling  wave  structures  appear  upstream  of  the 
cascade  inflow  hinting  at  the  existence  of  an  instability  mechanism  (Fig.  29).  As  for  the  turbulent  endwall 
boundary  layer  (Fig.  10)  the  endwall  region  near  the  pressure  side  of  the  blade  is  relatively  “quiet”.  Neither 
the  skin- friction  lines  nor  the  streamlines  hint  at  the  existence  of  a  pronounced  passage  vortex  (Fig.  30).  This 
conclusion  is  supported  by  visualizations  of  the  time-averaged  flow  (Fig.  31).  A  wall-tangential  streamwise 
vortex  that  is  about  one  third  of  the  blade  spacing  displaced  from  the  suction  side  of  the  blade  is  seen  in 
the  time-average.  A  similar  structure  can  also  be  observed  for  the  turbulent  endwall  boundary  layer  albeit 
with  much  reduced  coherence  (Fig.  12).  For  the  laminar  endwall  boundary  layer  the  vortex  originates  from 
a  region  where  the  skin-friction  lines  form  a  saddle  point  (Figs.  30  &  31). 
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Figure  29.  Iso-surfaces  of  Q  —  10. 


Figure  30.  Instantaneous  iso-surfaces  of  Q  =  10.  Red  lines:  Wall  skin-friction  lines  for  time  average.  Green 
lines:  Streamlines  at  z  =  0.05  for  time-average. 


A  comparison  of  Figs.  15  (turbulent  endwall  boundary  layer)  and  32  (laminar  end  wall  boundary  layer) 
reveals  that  for  the  latter  the  suction  side  trailing  edge  separation  is  reduced  and  the  pressure  side  separation 
is  almost  absent. 
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Figure  31.  Time-average.  Iso-surfaces  of  Q  =  10. 
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Figure  32.  Wall  pressure  coefficient  iso-contours  and  skin-friction  lines. 


Finally,  a  comparison  of  Figs.  16  &  33  and  Figs.  IT  &  34  provides  support  for  the  statement  that  the 
passage  vortex  is  virtually  missing  for  the  laminar  endwall  boundary  layer. 


Figure  33.  Time-average.  Iso-contours  of  vorticity  magnitude,  uout‘ 
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Figure  34.  Vorticity  iso-contours,  u)xt . 


For  completeness,  vorticity  contours  for  the  outflow  plane  are  shown  in  Fig.  35. 
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Figure  35.  Iso-contours  of  uo0ut  at  x=1.5Cx. 


2.  Blade  With  Fillet 

Finally,  the  L2F-EF  geometry  with  laminar  endwall  boundary  layer  was  investigated.  For  this  simulation 
the  fifth-order-accurate  discretization  for  the  convective  terms  was  chosen.  Instantaneous  iso-surfaces  of  the 
Q  vortex  identification  criterion  in  Fig.  36  &  37  reveal  a  hairpin  structure  upstream  of  the  blade  and  several 
similar  structures  downstream  of  the  leading  edge.  Also,  similar  to  the  simulation  for  the  turbulent  endwall 
boundary  layer  (Fig.  19),  streamwise  structures  can  be  observed  on  the  pressure  side  near  the  endwall. 
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Figure  37.  Instantaneous  iso-surfaces  of  Q  =  10.  Red  lines:  Wall  skin-friction  lines  for  time  average.  Green 
lines:  Streamlines  at  z  =  0.05  for  time-average. 


With  fillet  a  passage  vortex  is  originating  at  the  leading  edge  of  the  blade  (Fig.  38).  The  vortex  remains 
close  to  the  pressure  side  of  the  blade  and  dissipates  shortly  downstream  of  the  trailing  edge.  Interestingly, 
for  the  straight  blade  with  laminar  endwall  (Fig.  31)  the  passage  vortex  is  missing.  This  is  contrast  to  the 
results  for  the  turbulent  endwall  boundary  layer.  One  must  conclude  that  the  endwall  boundary  layer  state 
has  a  large  impact  on  the  mean  flow  field.  However,  caution  is  in  order  as  the  two  simulations  were  not 
carried  out  with  the  same  discretization. 
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Figure  38.  Time-average.  Iso-surfaces  of  Q  =  10. 


Figure  39.  Wall  pressure  coefficient  iso-contours  and  skin-friction  lines. 
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The  skin-friction  lines  and  pressure  iso-contours  in  Fig.  39  are  qualitatively  similar  to  the  results  for  the 
turbulent  endwall  boundary  layer  (Fig.  23).  The  suction  side  trailing  edge  corner  separation  is  narrower  for 
the  laminar  endwall  boundary  layer. 


Figure  40.  Time-average.  Iso-contours  of  vorticity  magnitude,  uo0ut> 


Iso-contours  of  the  vorticity  in  the  cascade  (Fig.  40)  provide  proof  of  a  strong  suction  side  corner  vortex. 
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Figure  41.  Laminar  endwall  boundary  layer  (grid  4).  Vorticity  iso-contours,  ujxi  . 


The  passage  vortex  is  visible  in  Fig.  41.  It  remains  close  to  the  pressure  side.  Iso-contours  of  the 
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Figure  42.  Iso-contours  of  ujout  at  x=1.5Cx- 


streamwise  vorticity,  ujout ,  at  the  outflow  plane  are  shown  in  Fig.  42. 
3.  Comparison 


Figure  43.  Skin-friction  lines,  wall-pressure  iso-contours,  and  streamlines  at  z=0.05.  Red  lines:  Turbulent 
endwall  boundary  layer.  Green  lines:  Laminar  endwall  boundary  layer. 

Direct  comparisons  of  the  results  for  the  straight  blade  with  turbulent  and  laminar  endwall  boundary 
layer  are  provided  in  Figs.  43  &  44.  The  boundary  layer  state  has  a  strong  effect  on  the  endwall  flow.  The 
most  prominent  differences  are  the  lack  of  the  passage  vortex,  a  reduced  suction  side  trailing  edge  corner 
separation,  and  a  more  two-dimensional  pressure  distribution  on  the  suction  side  for  the  laminar  endwall 
boundary  layer. 
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Figure  44.  Suction  side  skin-friction  lines  and  wall-pressure  iso-contours.  Red  lines:  Turbulent  endwall 
boundary  layer.  Green  lines:  Laminar  endwall  boundary  layer. 


C.  Total  Pressure  Loss  Coefficient 

Table  4  lists  the  mass- averaged  total  pressure,  dynamic  pressure,  and  total  pressure  loss  coefficient  for  the 
different  cases  as  well  as  for  the  RANS  calculations  by  Lyall  et  al.24  and  the  experiments  by  Sangston  et 
al.25  The  total  pressure  loss  coefficient  for  the  turbulent  endwall  boundary  layer  is  29%  (straight  blade)  and 
25%  (blade  with  fillet)  higher  than  for  the  laminar  endwall  boundary  layer.  This  once  again  illustrates  the 
strong  dependence  of  the  flow  on  the  endwall  boundary  layer  state.  For  the  turbulent  endwall  boundary  layer 
the  fillet  reduces  the  loss  coefficient  by  8.0%  (9.9%  reduction  in  experiment).  In  general,  the  loss  coefficient 
for  the  RANS  calculations  by  Lyall  et  al.24  and  the  experiments  by  Sangston  et  al.25  is  higher  than  for 
the  present  simulations.  This  may  be  explained  by  differences  in  the  postprocessing  of  the  data  (i.e.,  the 
computation  of  the  total  pressure  loss  coefficient),  by  a  lack  of  grid  resolution,  or  by  differences  w.r.t.  the 
endwall  boundary  layer  state. 


Case 

r" 

»,0,m 

n" 

V in 

r" 

u p,0,out 

Qout 

Ytot 

Straight  blade 

Turbulent  endwall  B.L. 

0.65876 

0.47128 

0.54411 

1.1148 

0.1147 

Laminar  endwall  B.L. 

0.47889 

0.46077 

0.38966 

1.0883 

0.08923 

RANS,  Lyall  et  al.24 

0.141 

Experiment,  Sangston  et  al.25 

0.142 

Blade  with  fillet 

Turbulent  endwall  B.L. 

0.65720 

0.47159 

0.55169 

1.1166 

0.1055 

Laminar  endwall  B.L. 

0.48802 

0.46246 

0.40381 

1.0908 

0.08421 

RANS,  Lyall  et  al.24 

0.131 

Experiment,  Sangston  et  al.25 

0.128 

Table  4.  Mass-averaged  total  pressure  coefficient,  c"  0,  and  dynamic  pressure,  q" .  Total  pressure  loss  coefficient, 

Ytot- 


IV.  Conclusions 

Three-dimensional  Implicit  Large  Eddy  Simulations  (ILES)  of  the  L2F  blade  with  and  without  fillet  were 
carried  out  for  a  Reynolds  number  based  on  axial  chord  of  Re=100,000.  Both  a  laminar  and  a  turbulent 
endwall  boundary  layer  were  considered.  The  results  for  the  straight  blade  (no  fillet)  with  turbulent  endwall 
boundary  layer  are  in  adequate  agreement  with  Particle  Image  Velocimetry  (PIV)  measurements  by  Sangston 
et  al.25  A  strong  passage  vortex  and  a  considerable  suction  side  trailing  edge  separation  are  observed.  Both 
are  total  pressure  loss  mechanisms  and  amenable  to  flow  control.  The  simulation  for  the  L2F  blade  with 
fillet  (L2F-EF  geometry)  and  turbulent  endwall  boundary  layer  indicates  a  greatly  reduced  coherence  of  the 
passage  vortex  in  accordance  with  the  experimental  findings.  Overall,  the  total  pressure  loss  coefficient  for 
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both  simulations  is  lower  than  in  the  experiment.  This  may  be  attributed  to  insufficient  grid  resolution. 
As  a  result  the  endwall  boundary  layer  is  too  thin.  Simulations  with  laminar  endwall  boundary  layer  (i.e., 
the  thinnest  possible  boundary  layer)  revealed  that  the  endwall  boundary  layer  has  a  profound  effect  of  the 
flow  topology.  For  the  straight  blade  (no  fillet)  and  with  laminar  endwall  boundary  layer  the  passage  vortex 
is  missing  and  the  total  pressure  loss  coefficient  is  greatly  reduced  (compared  to  the  case  with  turbulent 
endwall  boundary  layer).  With  fillet  a  weak  passage  vortex  is  obtained.  The  results  for  the  turbulent  and 
laminar  endwall  boundary  layer  are  diametrically  opposed:  For  the  turbulent  endwall  boundary  layer  the 
fillet  suppresses  the  passage  vortex.  For  the  laminar  endwall  boundary  layer  the  fillet  generates  the  passage 
vortex.  This  strong  dependence  of  the  flow  topology  on  the  endwall  boundary  layer  state  is  one  of  the  most 
interesting  outcomes  of  the  present  investigation. 

Overall,  the  mean  flow  topology  and  the  mechanisms  resulting  in  the  observed  unsteady  coherent  struc¬ 
tures  remain  poorly  understood.  One  of  the  objectives  of  the  project  was  to  employ  unsteady  Active  Flow 
Control  (AFC)  near  the  endwall  to  probe  instabilities.  Due  to  the  unexpected  complexity  of  the  endwall  flow 
simulations  and  the  simulation  turnaround  times  (a  typical  simulation  required  four  days  on  256  processors) 
the  planned  AFC  simulations  could  not  be  accomplished.  The  following  tasks  are  planned  for  the  future: 
The  dependence  of  the  results  on  the  grid  resolution  and  the  accuracy  of  the  numerical  discretization  will 
be  investigated.  The  computational  domain  for  generating  the  turbulent  endwall  boundary  layer  was  quite 
narrow  and  resulted  in  a  pitchwise  periodicity  of  the  endwall  boundary  layer.  The  narrow  approach  flow 
domain  will  be  replaced  by  a  wider  domain.  Simulations  at  a  reduced  Reynolds  number  (e.g.,  Re=25,000)  are 
less  costly  and,  should  the  physics  be  similar,  will  be  considered  for  investigating  the  endwall  flow  physics. 
Hybrid  turbulence  models  are  are  tempting  alternative  but  the  transitional  nature  of  the  flow  raises  questions 
with  respect  to  the  applicability  of  the  approach.  The  time-dependent  flow  data  will  be  analyzed  in  more 
detail  with  the  Proper  Orthogonal  Decomposition  (POD),33,34  and  based  on  the  analysis  different  Active 
Flow  Control  (AFC)  strategies  will  be  proposed. 
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